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The large-scale power deficit in the cosmic microwave background fluctuations might be rele¬ 
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inflation, which can provide a singular-free realization of inflation. We investigate the primordial 
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with different background dynamics, and give a universal formula for the power spectrum of pri¬ 
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to corresponding cases and show how the intensity of large-scale power suppression is affected by 
the pre-inflationary physics. 
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I. INTRODUCTION 


As the paradigm of the early universe, inflation has been generally regarded as a possible 
solution of the horizon, flatness, entropy, homogeneity, isotropy and primordial monopole 
problems [1],[2],[3],[4]. But maybe far more attractive is that inflation can generate the 
primordial perturbations, which have grown into all the structures observed in our universe 
today. The observations of cosmic microwave background (CMB) by Planck and WMAP 
has provided us with more and more information of the early universe, which shows that 
the single field slow-roll inflationary model is more likely to be the right one. 

However, a large-scale (or low-/) power deficit in CMB TT-mode spectrum observed by 
WMAP [5] and recently confirmed by the Planck Collaboration [6] [7] with higher precision 
is not concordant with the standard slow-roll inflation. It is hard to attribute this power 
deficit to the foreground as it has been observed by experiments with higher and higher 
statistical significance. Though the cosmic variance could be a source of this deficit, it is 
still very likely that the large-scale anomalies are induced by the physics preceding inflation, 
as the larger are the scales of the perturbations, the earlier are the times corresponding to 
their horizons exiting. 

The slow-roll inflation might last for just the minimal number of e-folds, i.e., just enough 
[8], and thus the power deficit on a large scale may be attributed to the pre-inflationary 
non-slow-roll evolution. In this case, the Planck best-fit single-held inhationary model only 
actually provides a ht for the intermediate and small angular scales. After the WMAP1 
data were released, some studies have been done in Refs. [9],[10],[11],[12],[13],[14],[15] along 
this line, and also recently [16],[17],[18],[19],[20],[21],[22],[23]; see Ref. [24] for a review. 

However, it might be more interesting that the large-scale anomalies could be relevant 
with the physics solving the initial singularity problem of inflation [11],[19]. In the bouncing 
model (see, e.g., [25],[26] for reviews), initially the universe is in a contracting phase, and 
then it bounces into an expanding phase, which results in a solution to the cosmological 
singularity problem. In Refs. [11],[14],[15],[16],[17],[21], it was noticed that if the universe is 
initially in a contracting phase and after the bounce it begins to inhate, the power spectrum 
of primordial perturbations will get a large-scale cutoff, which may naturally explain the 
power dehcit of the CMB TT-mode spectrum at low-/; see, e.g., [16] for details. In addition, 
it was also observed in [27] that if the pre-inflationary bounce actually occurs, the BB- 
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mode correlation at low-/ is also suppressed, while the TB- and EB-mode correlations on a 
corresponding scale may be enhanced. 

The superinflationary phase before slow-roll inflation may also provide a singular-free 
realization of inflation, which in the meantime explains the anomalies of the CMB power 
spectrum [18],[28]; see also [29],[30] for an almost flat pre-inflationary universe. 

In the above pre-inflationary scenarios, e.g., [11],[16],[18], initially the primordial pertur¬ 
bation is deep inside the horizon, which naturally set itself in Bunch-Davies(BD) vacuum. 
This implies that we can calculate a large-scale power spectrum without any assumption for 
the initial state of primordial perturbations. However, it is possible that the pre-inflationary 
era might consist of multiple phases with different background evolution, e.g., [24], In a 
certain sense, the introduction of multiple pre-inflationary phases may better simulate the 
physics of the pre-inflationary era, since due to the complexity of pre-inflationary physics, 
sometimes a single phase can hardly reflect the drastic change of the background parameters, 
e.g., [20]. Therefore, it is interesting to have a quantitative estimate for the power spectra 
of primordial perturbations from an arbitrary pre-inflationary era, involving the multiple 
phases with different background dynamics. 

This paper is organized as follows. In Sec.IIA, we introduce the evolution of the pre- 
inflationary background, which we will focus on, consisting of multiple phases. We require 
that the primordial perturbations can be produced in these phases. In Sec.IIB, we perform 
a model-independent calculation for the primordial perturbations and give a universal for¬ 
mula for the power spectra of primordial perturbations in terms of the recursive Bogoliubov 
coefficients. In Sec. Ill, we apply our formula to the bounce inflation and the superinflation 
preceding slow-roll inflation and show how the intensity of the large-scale power suppression 
of a primordial spectrum is affected by the pre-inflationary physics. We will see that a large- 
scale suppressed primordial spectrum may result in the power deficit at low-/ in the CMB 
TT-mode spectrum; however, the intensity of the power suppression is model-dependent. 
Sec. IV is the discussion. 

In Appendix A, we will investigate the Wronskian constraint for primordial perturbations 
and argue that in a certain sense the scenario with pre-inflationary era is equivalent to the 
inflation scenario with the non-Bunch-Davis initial state. In Appendix B, we will approxi¬ 
mately estimate the spectra index of primordial perturbations produced in each phase. 
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II. PRE-INFLATION ARY PRIMORDIAL PERTURBATIONS 


A. The pre-inflationary background 


The pre-inflationary era may consist of multiple phases. We define the inflation as phase 
0, the latest pre-inflationary phase as phase 1, and so on. 

The cosmological evolution of phase i is 

2 2 

or (—r]) 1+3ui i : (1) 


where u>i = Pi/Pi is a constant, u>i —1/3, and rj — f dt/a is the conformal time. The phase 
i may be expanding or contracting. 

However, noting the continuities of a and \H\, the evolutions of the inflation and phase i 
can be rewritten as 


ai(v) 


a o(v) = 


ao(vo 




1 + 


1 - Uinf ■ (V ~ T/o) ’ 

1 + • (■ rj - rji-i) 


l+3ct)j 


( 2 ) 

( 3 ) 


respectively, where rj 0 is that at the onset of inflation, = V-oivo) is the conformal 

Hubble parameter during slow-roll inflation, and r/,_ i is the matching time between phase i 
and phase i — 1, which signals the onset of phase i — 1. 

It seems that cu, is arbitrary. However, if we require that initially all perturbation modes 
are in the BD state, which is right only if its wavelength A 1/H, uii will be constrained. 
That the perturbation mode with A ~ 1/k extends outside the horizon in phase i marks 
the primordial perturbation that is produced in this phase, which requires %;(??*) < k < 
1 ). Thus \'Hi(r))\ must increase with time, which implies 


u>i < —1/3 for the expanding phase, 
or Ui > —1/3 for the contracting phase. (4) 


Or equally in other words, (4) must be satisfied to guarantee the primordial perturbations 
that can be produced during phase i. The slow-roll inflation is the evolution with oj ~ — 1, 
which satisfies (4). 

In general, there may be multiple pre-inflationary phases. Any one of them might be a 
contracting phase or an expanding phase. Therefore, there are two kinds of typical scenar- 
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ios which have aroused lots of interests, i.e., the bouncing scenario and the superinflation 
scenario, which we will focus on in Secs. III.A and III.B. 

In the bouncing scenario, a bounce happened near the end of the contracting phase then 
was followed by an expanding phase with u < —1/3. The bounce may be implemented 
by applying a higher-order derivative field [31],[32],[33],[34], which is ghost-free, and also 
viscous fluid [35] and modified gravity [36],[37],[38],[39],[40]. 

The expansion with oj < —1 is called superinflation, which is similar to the emergent 
scenario and describes a monotonically expanding universe with increasing energy density. 
The primordial perturbations generated during the superinflation have been studied earlier 
in Refs. [41],[42]. The case with co <C — 1 corresponds to the slow expansion scenario, which 
has been proposed in Ref. [44] (see also [45],[46] for the genesis scenario) and investigated 
in detail in Refs. [47],[48]. It is shown in Ref. [48] that there is no ghost instability during 
superinflation. Thus the pre-inflationary universe may also be superinflating. It is interesting 
to notice that both bounce and superinflation preceding slow-roll inflation may provide a 
singular-free realization of inflation. 

However, if the pre-inflationary era is the expanding phase with u — 1 characterized 
by fast-rolling dominance, e.g., [9],[49],[50],[51], except for the initial fast-roll inflation, e.g., 
[52], or the expanding phase with radiation dominance [10],[13], initially the perturbation 
mode should be outside the horizon, so it is not clear how to set its initial condition. It is 
possible that its initial value is set in a phase preceding fast-rolling dominance; however, this 
phase is still required to satisfy (4), e.g., a higher energy inflation preceding the fast-rolling 
phase [53]; see [24] for comments, and see also [23] for a domain wall or cosmic string phase. 
Here, we will not involve this issue and will assume that all pre-inflationary phases satisfy 
(4), which assures initially all perturbation modes are naturally set in BD vacuum and each 
phase may contribute to the production of primordial perturbations with the certain range 
of wave numbers. 

We qualitatively plot the evolutions of background and perturbation modes, which we 
will focus on, in Fig.l. 
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FIG. 1: The green line denotes the evolution of the comoving Hubble radius and r/i signals 
the onset of phase i. The blue dashed lines denote the evolutions of the primordial perturbations. 


B. The power spectrum of pre-inflationary perturbations 


The equation of the primordial perturbation 1Z is [54],[55] 


u" + 


cik 2 



( 5 ) 


where u = zTZ , the prime is the derivative with respect to conformal time, z = a-sjTM^kjc s 
[56], and z — a^2M 2 \e\/c s [48] for the super inflat ionary phase, and the definition of e is 
e = —H/H 2 = . We assume c 2 — 1 for simplicity. 

In the slow-roll inflationary phase, we have 


7 " 

z o 


Wl f 


z o [1 - T~L inf (t) - 770 )] 


2 > 


( 6 ) 


the solution of Eq.(5) is 


Uq 


~ V k^effo Co,l-^3/2 ( kVeffo) T @0, 2H3/2 ( &Veffo ) 


( 2 ) 


( 7 ) 


where and H^ 2 are the (3/2)th order Hankel functions of the first and second kinds, 
respectively, and r) e ff 0 = rj - rj 0 - 
In pre-inflationary phase i, we have 


£? __ (1 - 3ui)nUVi-i) _ 

* - 2 [1 + i±§^ni(vi-i) • (v - ^-i)] 2 ’ 


( 8 ) 
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the solution of Eq.(5) is 


Ui = 



( 9 ) 


where V{ 


3 1 -uij 

2 l+3aji 


and are the u*th order Hankel functions of the first and 


second kinds, respectively, and r/ e //. = rj — rji-i + (i +3cj .)^.(,^. ^ I 11 fact; Eq. (7) can be 
obtained from Eq.(9) while i — 0, but r]i_i should be replaced by r] 0 . 

Here, we require that around and at the matching surface besides there is not the ghost 
instability; there is also not gradient instability, i.e., (? s > 0. In this case, the perturbation 
can continuously pass through the matching surface between two adjacent phases, and its 
spectrum is insensitive to the physical details around the matching surface; see e.g., [57] for 
the bounce. The coefficients C % in Eqs.(7) and (9) are determined by requiring the continuity 
of u and v! at the matching surface. We can write the coefficients CVi and 2 of phase i 
recursively as 



M {i ’ i+1) x x • • • x M ( 



( 10 ) 


where the earliest phase is defined as the phase i ma x, and 


Mu M\2 
M 21 M 22 
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recursive matrix, which is given by 

iny/xy f 


Mu = 


\ \-H% vw (y) + J H^(x) + [ifhCM - »!?,,»] <>,(») 


+ (x '-y ')hW(x)hV> (y) f 


Mi 2 ~ 


m^/xy f 

8 l 


-H% v ty) + H™ ,(y) H S\*) + H^ +K (x) - H^(x) HgUy) 


r( 2) 


( 2 ) 


r( 2 ) 




^21 — 


iny/xy f 


H-L,Jy) - <‘’M - \ H -AW - x£l,(v) 


\ L 

(x~' - y~ 1 )Hl]\x)Hi,]l I (y) 1 


A^22 — 


m^Jxy f 


, //( 2 ) 


M-<’„„(») <’(*)- - H^x) H^Jy) 


(1) 


r(l) 


(x x - y l )H£\x)HW +i (y)\, 


(11) 

where x = — 2 (y. , , and y = —% 7 —n- When i = 0, r+i should be replaced 

(1+3 Ldi)Hi(rn) (l+3wi + i/Hi + i(?7i) ’ 11 1 ^ 

with rjo- A result similar to (11) was obtained in [24], 

In the earliest phase, the coefficients Cl i and Cl 2 are determined by the initial 

JL 7 ('■maxi 1 - L maxi^ J 

z" 

condition. When k 2 lmax , the perturbation mode is deep inside the horizon, which is set 

^i-max 

in BD vacuum, 


u ~ 


V2k 


0 —ikr) 


( 12 ) 


When k 2 lrnax , Ui max given in Eq.(9) should approximate to the form in Eq.(12). Thus, 

Zimax 

we get 


a 




2v^ 


5 ^^max j2 


(13) 


Obviously, Cl mtw;) i and Ci maX} 2 satisfy the so-called Wronskian (or canonical normalization) 
constraint (see e.g. [58], [59]) 


4 k 


— (Id i| 2 -|Cl 2 | 2 )=i. 

7 r 


(14) 


And actually, for the following phase i, C)^ and Cl^ will always satisfy the Wronskian 
constraint, which will be proved in Appendix A. C h \ and Cl^ are related to the so-called 
Bogoliubov coefficients by Eq.(25). 

The power spectrum of 1Z is 


'Pn — 


k 3 

2^ 


Mo 

*0 


(15) 
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Hi n f 


•C 1, we get a uni- 


After substituting Eq.(7) into Eq.(15) and requiring 
versal formula 


V n = V% f -k |C 0) i -C 0 , 2 \ 2 , (16) 

71 

where V^ f = 2M i e . f is the spectrum of the standard slow-roll inflation, H m f = 

H 0 (vo)/a 0 (vo) is the Hubble parameter during inflation, and in this sense actually 'Ho(bo) = 
'Hinf. The spectral index of 1Z is 


^7 Z ^inf “I - 


dln(fc|Co, 1 -C , o, 2 | 2 ) 

dink 


( 17 ) 


where rii n f = —yjpf—f 1 is the spectral index of slow-roll inflation, which is nearly unity. 

The coefficients Cop and Co, 2 are determined by the recursive Eq.(10), and thus the 
effects of all pre-inflationary phases are encoded in Cop and Co )2 , which are nontrivial. By 
calculating Eq.(16), it can be found that the perturbations produced in phase i roughly have 
a power spectrum 

V n ~k 3 ~ 2vi . (18) 


We give a proof for this result in Appendix B. It is observed in Appendix B that the power 
spectrum of perturbations is Eq.(18) but modulated with a small oscillation, which is induced 
by the evolution of perturbation through the matching surface between adjacent phases. 


III. THE LARGE-SCALE POWER SUPPRESSION 


In this section, we will apply our universal formula (16) to the bounce inflation and the 
superinflation preceding slow-roll inflation, and show how the intensity of the large-scale 
power suppression in the CMB fluctuations is affected by the pre-inflationary evolution. 

The large-scale power suppression requires that the e-folding number of slow-roll inflation 
is just enough or less, and in the meantime the pre-inflationary era can contribute a strong 
blue-tilt spectrum. The spectral index of pre-inflationary perturbations is approximately 




dln(fc|C 0 ,i-C 0 , 2 | 

d. In k 


(19) 


Thus the intensity of suppression can be model dependent, since |C 0 ,i — (Tog1 2 is different 
for different parameters cuj and r\i. 
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Here, we will focus on the cases that the pre-inflationary era consists of two phases, i.e., 
imax = 2 , and set Poivo) = Pmf and r] 0 = 0 is the time at the onset of inflation; thus, we 
have 77 < 0 in the pre-inflationary phases and 77 > 0 in the inflationary phase. According to 
Sec.IIB, for i max = 2, the solutions of Eq.(5) are 

u 0 = 

Ml = 


U 2 = 

respectively, where 771 < 0 is the conformal time at the matching surface between phase 1 
and phase 2 . 

We define 

1i,ui,uj 2 ) = — 777 j = — k\C 0 i — C 02 I 2 , (21) 

7T 

where according to Eq.(10), we have 



(!) / I_, k \ , n w( 2 ) t ’ k 


Co,iH 3 / 2 (—kr) + ——) + C 0 , 2 H 3 ^(—kr) + 

tl'inf Jl'inf 

2k 


(i + 3u 1 )n 1 (o)'- 


C^HiH-krj- 


+C lfl H^(-kn- 


(1 + 3^ 1 )^ 1 (0)' 
2k 


(1+ 3070^(0) ■ 


7 r 


2 v V + Vl (l + 3u 2 )n 2 (vi^ V2 


H^X-krj + /ct7i - 


2k 


(1 + 3u 2 )'H 2 (r] 1 


-), ( 20 ) 


Q 

C, 


^ = M ( °P x x 


0,2 



( 22 ) 


and the components of A4^ 0,1 \ A4^ 1,2 ^ can be obtained from Eq.(ll). P(r] 1 , u 1 , u: 2 ) reflects 
the shape of the spectrum. We will check how the shape of the spectrum, as well as the 
intensity of suppression on a large scale, changes with the parameters 771 , uj\, and uj 2 in the 
bounce inflation and the superinflation preceding slow-roll inflation. 


A. The superinflationary phase before inflation 

The superinflation is defined as the evolution with e = —H/H 2 < 0 or oj < —1, e.g., 
[41],[42],[43],[44]. The model with one single superinflationary phase oj = —5/3 preceding 
slow-roll inflation has been built in string theory in Ref. [18]. In Refs. [29],[30], the pre- 
inflationary universe is in a slowly expanding genesis phase with e -C — 1 , which is almost 
Minkowski space. This genesis phase actually also belongs to the superinflation, but since 
e <C — 1 , the expansion is actually slow; see also [28] for the emergent universe. 
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We will fix phase 2 with u >2 = —5/3. Thus in a certain sense phase 1 actually corresponds 
to an intermediate phase, which obviously must also be expanding. The duration that this 
intermediate phase lasts is A rj = l^l- Here, the parameter uq in phase 1 is model-dependent 
but should satisfy uq < —1/3. By requiring the continuity of \H\, we get 


ftl(O) = 'Hinf , 




_ 2 T-Ljnf _ 

2 + (1 + 3ui)rjiHi n f 


(23) 


where Hmf is determined by the amplitude of CMB fluctuations. Thus in equation set (20) 
only parameters rji and uq are left to be free. 

We plot P(rji, oq, — 5/3) in Fig.2 with different iji and uq. The black solid curve, i.e., 
P(0, —1, —5/3), is the case with only one single superinflationary phase before slow-roll 
inflation, which has been studied in Ref. [18]. The perturbation mode with wavelength 
1/k — 1 /'Hmf is that exiting the horizon at conformal time rj = 0, i.e., the onset of inflation. 
In Fig. 3, we plot -P(r/i,uq, —5/3) by replacing k with k/\rji\, which corresponds to move the 
suppression of the spectrum to a smaller scale. 

In Fig.2, the effects of an intermediate phase between superinflation and slow-roll inflation 
is obvious, compared with the case without the intermediate phase. The first peaks in the 
dashed curves will left shift with the increase of the duration \r]i\ that the intermediate 
phase lasts. The height of the first peak will lower with the decrease of uq, since the first 
peak corresponds to tt— < 1 while the spectrum is blue tilted. The power spectrum of 
perturbation from the intermediate phase , i.e. phase 1, is roughly Vn ~ (tt—) 3 ~ 2vi , which 

v '•'inf 

will be proved in Appendix B. 


B. The bounce inflation 

The pre-inflationary universe may be contracting. After the bounce, the slow-roll inflation 
begins, which is the so-called bounce inflation scenario [11],[14],[15],[16],[17], [21] with G- 
bounce, and [22] with the quintom. It has been found that in this scenario the power 
spectrum of primordial perturbations will get a large-scale cutoff, which may lead to the 
power deficit of CMB TT-mode on a large angular scale. 

However, around the bounce the evolution of the universe is complicated. It might 
be possible that after the bounce the universe enters into an intermediate phase prior to 
the slow-roll inflation. We will check how the shape of the spectrum changes with this 
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FIG. 2: The power spectrum P(r/ 1 , wi, —5/3) of the model, in which there is an intermediate phase 
between superinflation and slow-roll inflation, with different ij i and oj\. P( 0, —1, —5/3) corresponds 
to the case in Ref. [18] without the intermediate phase. 








FIG. 3: P is obtained by replacing k with k/\r/i\. 


intermediate phase. We fixed U 2 = 1 and assume that the intermediate phase is characterized 
by the constant state parameter cci, which is model dependent but satisfies uj\ < —1/3. By 
requiring the continuity of |"H|, we get 


ftl(0) = Hinf , 


^(Vl) 


27-Linf 


2 + (1 + 3ui)r]i'H in f 


(24) 


where is determined by the amplitude of CMB fluctuations. Thus in equation set (20) 
only parameters 771 and oj\ are left to be free. 

We plot P(/ 7 i,o;i,l) with different 77 ! and u 1 in Fig. 4, in which the black solid curves, 
i.e. P(0, —1,1), correspond to the bounce inflation studied explicitly in [11] [16]. In Fig. 4, 
the effects of an intermediate phase between contraction and slow-roll inflation is obvious. 
The case is similar to IIIA. In Fig.5, we plot P(r]i,uji, 1) by replacing k with k/\r)i\, which 
corresponds to move the suppression of the spectrum to a smaller scale. Actually, similar 
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FIG. 4: The power spectrum P(t)i,uji, 1) of the model, in which there is an intermediate phase 
in the bounce inflation scenario, with different r]\ and u)\. P(0, — 1,1) corresponds to the case in 
Ref. [18] where there is only one single contracting phase before inflation. 




FIG. 5: P is obtained by replacing k with k/\rji\. 
shapes of the power spectrum have also been obtained in Ref. [20]. 


IV. DISCUSSION 

The large-scale power deficit in the CMB TT-mode spectrum may imply certain pre¬ 
inflationary physics, e.g., a contracting phase followed by the bounce or a superinflationary 
phase before slow-roll inflation, which can provide a singular-free realization of inflation. 
However, the physics of the pre-inflationary era might be complex; sometimes a single phase 
can hardly reflect the drastic change of the background dynamics. Thus it is interesting to 
have a quantitative estimate for the power spectrum of primordial perturbations from an 
arbitrary pre-inflationary era, involving multiple phases with different background dynamics. 
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We perform a model-independent calculation for the power spectrum of primordial pertur¬ 
bations produced during the pre-inflation with different background evolutions. We require 
the relevant physical quantities to continuously pass through the matching surface between 
adjacent phases, and we obtain a universal formula (16) for the primordial spectrum in terms 
of the recursive Bogoliubov coefficients. 

We apply our formula to the bounce inflation and the superinflation preceding slow- 
roll inflation and show how the intensity of the CMB power suppression on the large scale 
is affected by the pre-inflationary physics. It is found that due to the existence of the 
intermediate phase, the intensity of the power suppression becomes model dependent. 

We, with the power spectrum in Fig.4, plot the CMB TT-mode spectrum in Fig.6, in 
which Planck 2013 data are used and the model with an intermediate phase is called the 
extended bounce inflation. It has been found in Ref. [16] that the bounce inflation model 
can improve the fit to the data with ps —4.6 with respect to the standard inflation 

model with the power-law spectrum. Thus, it is interesting to have a global fitting analysis 
with the Planck new data, which is left in the upcoming work. 

In addition, both the bounce inflation and the superinflation preceding slow-roll inflation 
may also explain a large dipole power asymmetry at low -l in CMB TT-mode spectrum 
[16],[18]; see also other attempts [60],[61],[62],[63],[64],[65], as well as [66],[67]. Actually, 
there may also be dipole power asymmetry in CMB polarization [68],[69], which might be 
larger than those in the TT-mode power spectrum. Thus it is interesting to have an estimate 
for the dipole power asymmetry of CMB polarization in the scenarios discussed here. 

For the sake of simplifying the analysis and algorithm, we have made several assumptions. 
We have divided the pre-inflationary evolution into multiple phases and assumed that each 
phase possesses a constant equation of state parameter, and the transitions from one phase 
to the next are instantaneous. But in reality, the equation of state parameter may be 
nonconstant, the transitions should be smooth, and the equation of state parameter should 
change smoothly. However, as long as there are long phases with approximately constant 
equation of state parameters and relatively quick transitions between the different phases, the 
approximation we adopted to simplify our analysis and algorithm is reasonable. Of course, 
the influence on the perturbations induced by a nonconstant equation of state parameter 
is interesting, and the effects of the modes that exit the Hubble radius near the transition 
from one phase to another are even more intriguing, which have been studied in [70],[71]. 
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FIG. 6: The CMB TT-mode power spectra for inflation with power-law power spectrum(red solid 
line), bounce inflation model(green dashed line), and extended bounce inflation(blue dashed line), 
in which Q b h 2 = 0.02203, Q c h 2 = 0.1204, r = 0.09, 100© s = 1.04, A, = 2.1 x 10” 9 , n s = 0.961. 
The orange points are the Planck 2013 data with la errors. 

Throughout, we have used perturbation and background equations that are valid for 
general relativity (GR), since we lack the knowledge of the new physics that is needed for 
the pre-inflationary evolution. To implement bounce and superinflation evolutions in GR, 
we must violate the null energy condition (NEC) or require a closed universe. Otherwise, 
one must go beyond GR. 

The violation of NEC in GR usually leads to ghosts, which indicate dangerous instability. 
One of the several ways to avoid ghosts is by introducing the higher-order derivative scalar 
field. In some ghost-free Galileon bounce or superinflation models, by delicately designing 
the Lagrangian, one can get the perturbation equation similar to ours in the form (see, e.g., 
[31],[48],[72]), in which the sound speed is constant, which is set as 1 here. In that case, 
our analysis is still able to approximately capture the bounce or superinflation dynamics. 
Additionally, in singular bounce models, e.g., the original ekpyrotic model, the background is 
GR and the perturbation equation does not need to be modified, until the scale factor a(t) is 
so small that quantum gravity effects become important. Then, if the boirnce period is short 
so that the link between the contracting phase and the expanding phase is approximately 
instantaneous just as we have assumed, our analysis is also approximately valid. However, in 
some kinds of modified gravity theories, in which the perturbation equations are modified, 


15 








especially when higher-order terms of gravity appear in the perturbation equations, e.g., 
[73], our analysis will no longer be robust and will need to be revisited. 

It is significant that recently in [74], some general results for the evolutions of perturbation 
modes during bounce have become available, which might be used to track the perturbations 
during the bounce. And obtaining the recursive matrix for some specific bouncing model to 
better understand how the modes evolve during this period is an interesting issue, which we 
will back to in the future. 

We have assumed c 2 — 1 after Eq.(5), which is suitable for the cosmology driven by the 
scalar held. However, for an ideal fluid, (? s = oj, and for u < 0, it is negative which indicates 
an instability in the system. It is well known in thermal physics that a negative heat 
capacity typically implies that one is looking at thermal fluctuations around an incorrect 
(possibly tachyonic) vacuum. For oj < 0, essentially the short wavelength fluctuations grow 
exponentially rather than oscillate. Actually, for scalar fields, c 2 may also be different from 
1. The effects of varying sound speed on primordial perturbations has been studied; see, 
e.g., [75],[76]. But as long as c 2 is a positive constant during each phase, our result is not 
qualitatively altered. 

When deriving (16), we assume that around and at matching surfaces the perturbation 
mode has no ghost and gradient instabilities. How to implement such a requirement is an 
interesting issue, e.g., [30],[57]. However, if this requirement is not satisfied, which actually 
is not allowed physically, the result of the perturbation spectrum will be strongly affected 
by the physical details around the matching surface. 
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Appendix A: The Wronskian constraint for Q i and 2 

Initially k <C H, and the perturbation should be in its minimal energy state, i.e., BD 
vacuum. Thus, the BD vacuum is usually regarded as the initial state of the primordial 
perturbations. The effects of non-BD initial states on the primordial perturbation has been 
discussed in, e.g., Refs., [77],[78],[79],[80],[81],[82],[83],[84], Here, the Bogoliubov coefficients 


16 


are 


CP 2 \j Ci'h 

k 


A = 2 \I-C it2 . 


( 25 ) 


The Wronskian (or canonical normalization) constraint requires |«j| 2 — |A| 2 = 1- The initial 
state, i.e., BD vacuum, corresponds to cti max = 1 and /3i max = 0, i.e., 


a 


7r 


imax 5 1 


2 y/k' 


a 


imax 


= 0. 


(26) 


We will prove that if Ci maXt i and Ci max: 2 satisfy the Wronskian constraint, for any phase 
i, and C \ t 2 also satisfy this constraint. This is equivalent to the statement that for any 
phase i, if 


Ak 


7T 


(iQ+iql 2 — | Ci+1,21 2 ) — 1 


is satisfied, we always have 


It equals the proof that 


Ak 


— (|C' i ,i| 2 -|C ii2 | 2 ) =1. 


7 r 


Ci, 1 Ci 2 ) X I 1 0 ) x ( ^ I = I C w C i+1> 2 


0 -1 


a 


i ,2 


(27) 


(28) 


XI 1 ° x | (29) 

0-1 \ 0 + 1,2 


where C is the complex conjugation of C. Because of 


we have 


^' l ' 1 ‘ _ 1 ) x • 


C ia 


Ci 


7+1,2 


Ci,i C12 ) = | Cl>1 ) = ( Cl+1,1 \ X (A^ (M+1) ) f 


C i>2 J y Ci + 1,2 

= (c w C w ) x (A 4 (i,i+1) ) f . 

Thus after substituting Eqs.(30) and (31) into Eq.(29), we have 


(30) 


(31) 


x [ 1 ° j x (Al ( i i+1 )) = I 1 ° 


(32) 
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where Ml^ 1 ' 1 ^ is defined in Eq.(ll). 

Therefore, it is left to prove 

' \Mll \ 2 - \M. 2 l \ 2 = 1 
Mil X Ml2 — M.21 x M.22 — 0 

— — • ( 33 ) 

Ml 12 X Mi ll — Ml 22 X Ml 21 = 0 

k l-^d^l 2 — \Ml-22 1 2 = —1 

This actually can be obtained by using 

#F«) = fff K), (34) 

and 

Hl'HOH i 2 i> + „«) - (35) 

ffence, for any phase i, if (27) is satisfied, we always have (28). 

The Wronskian constraint for C' 0 j and Co ,2 is 

— {\C 0 ,i\ 2 -\C q , 2 \ 2 ) = 1. (36) 

7r 

Thus, Viz in (16) becomes 

Viz = Vn f h ICV - e 0 , 2 | 2 = V% f |a 0 - VI 2 , (37) 

where /3o 7 ^ 0. Thus, in a certain sense the scenario with the pre-inflationary era is equivalent 
to the inflation scenario with the non-BD initial state. 


Appendix B: The perturbation spectrum from pre-inflationary phase i 


In this appendix, we will prove that the power spectrum of perturbations produced in 
phase i is approximately 


Viz 


k 


n 


inf 


ni —1 


where 


and Vi is defined as v t = | 


Hi - 1 = 3- 2 v u 


(38) 


(39) 


1 —Uj 


l + 3tdj 
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First, we give a proof for the case with i max = 1. The wave number of the perturbations 
produced in phase i satisfies k <C 'K% n f. We assume that phase 1 is an expanding phase, 
thus oji < —1/3. According to Sec. II, we have 


u o 


u 1 


+ XT— (C o,i Hy 2 (-kr] + ——) + C Q ^HfJ 2 {—kr] + xt—) 
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—krj — 


2k 


(1 + 3u;i )'Hi n f 




2k 
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(1 + 3CU1 )'Hi n f 
2k 

(1 + 3w 1 )^ in/ V’ 


(40) 


where 

Cl,1 = 2 ^F c ' ll2 = 0 - (41) 

Because always appears with k as k/Hmf below, we will set = 1 for convenience, 
and actually can easily get it back. We have 


_ qjinf 
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(42) 


where 
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' 1 + 3cui 


1,1 v l + 3cV + 2kHl+vlK l + 3u 1 ' 


Hi, 


The special case is oj\ = —5/3, which gives V\ — 1, — 1 + Vi — 0. When k <C 1, i.e. 
<C 1, we get 


' n ' n nU 
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If —5/3 < cui < —1/3, then v\ > 1, — 1 + v\ > 0. When t r— <C 1, we get 

' ' 't-in f 
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( 46 ) 


If cui < —5/3, then v\ < 1, — 1 + v\ < 0, H^\ +Vl should be replaced by e*( 
When ^r— C 1, we get 
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Thus, we always have Vr 
( jj\ > —1/3. 


Hi. 


Tli — 1 


. This result can also apply to the case with 


Then, we focus on the case with i ma x > 1- Though the spectra index seems nontrivial, 
based on some assumptions, we could get a similar result. 

We assume that the pre-inflationary era is phase i dominated, i.e., the phase i lasts long 
enough, nearly, (t^_ 1 — rf-) —> oo. The modes of the perturbations produced in phase i satisfy 
'H~ i (ry) <C k <C Thus, for i > i, we always have x 1 and y 1, x and y are 

those in Eq.(10), and we obtain 
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AI 2 I = 

.M 22 = 


x 1 — y 1 
4 i 


,-i(x-y — 


7) ~ e -i(x- V -2-£ 


T ) 


xX ~y 1 c -i( X+ y^J±H± l v) 


0 , 


xl -y 1 c i(x+y- v -^,) n 0) 


x 1 — y 1 


1 + - | e i(x-y-^p±*) _ e i(x-y-2-2±L*) 

4 i 


(48) 


20 


























Thus, when i > i, the recursive matrixes .A+®’® +1 ) are diagonal. We get 
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where E is real, C irnax: i = ^=. Because we are concerned about |<^ 0,1 — C+l in the final 
result, e~ l ' E makes no difference, we can ignore it and write that 


Q,i 

a ,) ^ \ 0 




(50) 


For i < i, we always have x -C 1, y <C 1. By expanding the Hankel functions to second 
order, we have 

Ci,l + Cj,2 = (-Mu + + (Ati2 + A / l22)Cj + i i 2 

/ fc \ Vi + 1 ~ Vi 

~ A,i' (Q+ip + 2 ) 
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where the constants A,i, A, 2 , A, 3 , A ,4 are independent with k, and Vi is defined as 
Because k <C < 'H{rj~ i _ 2 ) < ■■■ < T-L(rj 0 ), we will not care about 


Vi = 


l —m 


1+3 uji 


1 ), ...,"H(?7o) below, and just regard -C 1 as k <C 1 instead. 

According to Eqs.(50) and (51), we get 
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and then 
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( 53 ) 



After neglecting the high order terms, we have 



(54) 


where no = 3/2 for the inflationary phase. Therefore 



(55) 


It should be noted that (55) is a good approximation only if phase i lasts long enough, i.e., 


<C k <C The reason is apparent, since if phase i lasts only a very short 

time, the oscillations of spectrum around the matching surface would destroy the relation 
we proofed above. 

We plot rii — 1 with respect to cu* in Fig. 7 (see also [85]), which provides guidance for 
building a pre-inflationary model leading to a cutoff spectrum on a large scale. We see 
that the power spectrum is strong blue-tilt only for the expansion with cu; < — 1 and the 
contraction with uii > 0. 


4 



FIG. 7: The figure of rii — 1 with respect to Wj, which is plotted based on n* — 1 = 3 — 2 Vi and 


i-j ■ — if — ^4 
— 2 l+3u>i 


The magenta curve is rii ~ 1- And the blue dashed line is uji = —1/3. 
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